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ABSTRACT 

We present the last chapter of a spectroscopy program aimed at deriving the redshift or a lower 
limit to the redshift of BL Lac objects using medium resolution spectroscopy. Here we report new 
spectra for 33 BL Lac object candidates obtained in 2008-2009 confirming the BL Lac nature of 25 
sources and for 5 objects we obtained new redshifts. These new observations are combined with 
our previous data in order to construct a homogeneous sample of ~ 70 BL Lacs with high quality 
spectroscopy. All these spectra can be accessed at the website http://www.oapd.inaf.it/zbllac/. 
The average spectrum, beaming properties of the full sample, discussion on intervening systems 
and future perspectives are addressed. 

Subject headings: BL Lacertae objects: general. Galaxies: quasars: absorption lines. Instrumentation: 
spectrographs 
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1. Introduction 



2. Observations and data reduction 



BL Lac objects are extragalactic sources for 
which the redshift z cannot be easily derived. 
This is due to the combination of intrinsically 
weak emission lines and strongly beamed nuclear 
non thermal emission. The determination of z 
is obviously important in order to properly char- 
acterize their physical properties. While other 
distance indicators were introduced like the pho- 
tom etry of the host galaxy ( imaginK redshift, see 
e.g. 



Sbarufatti et al. 



(l2005aD'). or TeV 7 - ray a b- 
sorption by the EBL fe.g. iPrandini et all ( 2QlCll )). 
the direct spectroscopic line observation remains 
the only way to produce a robust measure of their 
distance. The determination of z is becoming 
even more compelling by the recent development 
of 7-ray astronomy, which demonstrated that the 
extragalactic 7-ray sky is dominated by this class 
of objects. 

In order to improve the quality (in terms of 
throughput, spectral resolution and S/N ratio) 
of the spectrum, we have carried out systematic 
campaigns at European Southern Observatory 
Very Large Telescope (VLT), aimed to observe 
~ 90 BL Lac candidates of unknown redshift. Re- 
sults on the fir s t three runs were presented in 
Sbarufatti et~al] (|2005bl E)06a. 2009 ). Paper I, II 
and III respectively , where the criteria of selec- 
tion of the sample (see also Section 15. ip , details 
on the used spectrometer, and the reduction and 
analysis procedure are given. 

In this paper we report on observations of 33 
new BL Lac candidates that complete our obser- 
vational campaign and for 25 of them we confirm 
their BL Lac nature. In Section [2] we summarize 
the observations and data reduction process and 
we show the results in Section [3l The average 
spectrum of the class is presented in Section 
while a discussions of the properties for the full 
sample is reported in Section [5l 

Throughout the paper we adopt the following 
concordant cosmology: Hq = 70 km s~^ Mpc~^, 
= 0.27 and 17a = 0.73. 



The spectra were obtained in Service Mode 
using the FOcal Red ucer and a low dis persion 
Spectrograph fF0RS2. ISeTfert et al.l tOO(t ). The 
observations took place between April 2008 and 
June 2008 at VLT-UTl (Antu) and from Septem- 
ber 2008 to February 2009 at VLT-UT2 (Kueyen). 
The spectrograph was configured with a 300V+I 
grism with a slit of 2" achieving the spectral range 
of 3800 - 8300 A with a dispersion of - 3. 2 A /pixel. 
Detailed informations on the observations of each 
source are given in Table [TJ The seeing during the 
observation was in the range 0.65"-2.50" with an 
average of 1.20" (rms - 0.4"). 

Data reduction was executed using iraeIiI adopt- 
ing standard procedures. In particular, for each 
target we performed bias subtraction, fiat field- 
ing and bad pixels or cosmic rays rejection. For 
each object we secured 2-6 individual exposures 
that were then combined into a single average 
image. The wavelength calibration was achieved 
using the spectra of a Helium Neon Argon lamp 
reaching an accuracy of ~ 3 A rms. Even though 
this program did not require photometric condi- 
tions during observations, checking the weather 
database at Paranal, the sky was clear during 
most of the time. This enables us to fiux cal ibrate 
our spectra using standard stars (|Okelll990l ). The 
fiux calibration was assessed through the relative 
photometric calibration of the acquisition image 
and the overall uncertainty is Am ^ 0.1 mag 
(rms). Finally, each spectrum has been dered- 
dened applying the extinction law described in 
Cardelli et al.l (11989) assum i ng th e Eb-v values 
computed by ISchlegel et al.l (|l998[ ) . 



3. Results 

From the spectroscopic point of view we con- 
sider as bona-fide BL Lac the sources that either 
have a featureless spectrum or in which their in- 
trinsic features have equivalent width (EW) < 5 
A. Based on this classification we confirm the 
BL Lac nature for 25 objects. Their spectra 



-'l IRAF (Image Reduction and Analysis Facility) is dis- 
tributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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Table 1 
Journal of the observations 





R.A.(J2000) 


Dccl.(J2000) 


R 


Date 


xi^xp osui c Tunc 


N 


S /N 




Rcf 






(2) 




(3) 


(A) 














GC 0109+224 


01 


12 09.103 


+ 22 


44 35.052 


15.60 


2008-09-20 


2040 


6 


200 


2.00 


vol 


PKS 0140-059 


01 


42 40.810 


- 05 


43 26.472 


18.20 


2008-10-26 


1362 


2 


65 


1.00 


V06 


RES 0337 


02 


38 29.549 


- 31 


16 33.312 


16.40 


2008-09-30 


1362 


2 


130 


0.65 


vol 


PKS 0302-623 


03 


03 56.678 


- 62 


11 40.812 


17.90 


2008-09-27 


1362 


2 


150 


1.40 


894 


RES 0499 


03 


58 55.903 


- 30 


54 02.664 


17.60 


2008-09-20 


1362 


2 


40 


1.00 


vol 


PKS 0410-519 


04 


11 41.858 


- 51 


49 31.548 


16.50 


2008-10-01 


1362 


2 


70 


1.60 


G04 


PKS 0439-299 


04 


41 22.987 


- 29 


52 31.296 


18.00 


2008-10-25 


1362 


2 


30 


1.20 


vol 


PKS 0443-387 


04 


45 14.474 


- 38 


38 34.296 


17.00 


2008-10-26 


1362 


2 


140 


1.70 


G04 


RES 0589 


04 


48 40.733 


- 16 


32 51.684 


16.60 


2008-10-28 


1362 


2 


130 


0.80 


vol 


ZS 0506-1-056 


05 


09 25.900 


-1- 05 


42 21.000 


15.10 


2008-11-02 


800 


2 


200 


2.00 


vol 


PMN J0529-3555 


05 


29 40.800 


~ 35 


55 41.160 


18.90 


2008-10-29 


1362 


2 


30 


1.30 


G04 


PKS 0539-543 


05 


40 50.200 


- 54 


17 58.700 


17.30 


2008-10-29 


1362 


2 


70 


1.60 


G04 


PKS 0607-605 


06 


08 00.334 


~ 60 


31 27.732 


17.60 


2008-10-29 


1362 


2 


50 


1.90 


G04 


MS-0622. 5-5256 


06 


23 42.667 


~ 52 


58 01.704 


19.20 


2008-10-26 


1362 


2 


45 


2.50 


vol 


PKS 0735-f 17 


07 


38 08.100 


+ 17 


43 03.600 


15.80 


2009-01-13 


2720 


4 


160 


1.30 


vol 


OJ-131 


08 


20 55.562 


- 12 


59 35.340 


16.50 


2008-04-09 


1362 


2 


100 


1.60 


vol 


RX J1022. 7-0112 


10 


22 46.704 


- 01 


13 11.280 


16.90 


2008-04-09 


1361 


2 


100 


1.50 


vol 


RX J1023. 9-4336 


10 


23 56.446 


- 43 


36 46.980 


16.70 


2008-04-10 


1361 


2 


200 


1.40 


vol 


PKS E1056-113 


10 


59 09.478 


- 11 


34 36.516 


16.20 


2008-04-18 


1361 


2 


100 


1.50 


V06 


lES 1106-1-244 


11 


09 18.245 


+ 24 


10 45.156 


18.70 


2008-04-22 


2720 


4 


40 


1.30 


vol 


RX J1117. 04-2014 


11 


17 07.522 


+ 20 


14 49.632 


15.10 


2009-02-16 


2612 


4 


130 


1.00 


vol 


RX J1241. 8-1455 


12 


41 49.862 


- 14 


56 42.756 


17.70 


2008-04-22 


2720 


4 


140 


0.95 


vol 


lES 1239-1-069 


12 


41 51.314 


+ 06 


36 06.840 


18.50 


2008-04-17/22 


2720 


4 


50 


1.00 


vol 


4C 12.46 


13 


09 31.862 


+ 11 


53 50.676 


18.20 


2008-04-18 


1361 


2 


70 


1.00 


vol 


OQ-240 


14 


27 03.307 


+ 23 


47 39.840 


14.50 


2008-04-07/17 


800 


4 


170 


1.00 


vol 


OR 103 


15 


04 27.670 


+ 10 


29 17.520 


18.90 


2009-02-16 


1361 


2 


150 


1.70 


vol 


RES 1478 


15 


16 16.718 


- 15 


24 18.000 


18.50 


2008-04-02 


1361 


2 


50 


0.90 


vol 


PKS 1815-553 


18 


19 44.678 


- 55 


22 06.204 


21.50 


2008-04-22 


1362 


2 


10 


1.30 


894 


OW-080 


20 


50 08.688 


+ 04 


08 15.792 


19.00 


2008-04-22 


1361 


2 


90 


0.70 


vol 


OX-183 


21 


52 27.972 


+ 17 


34 30.396 


18.00 


2008-06-02 


2722 


4 


50 


0.90 


vol 


RES 1899 


22 


49 13.829 


~ 13 


00 05.616 


18.20 


2008-06-08 


1361 


2 


60 


1.00 


vol 


RX J2319. 6-1-1611 


23 


19 44.371 


+ 16 


12 32.688 


17.30 


2008-06-08 


1361 


2 


100 


1.70 


vol 


lES 2322-409 


23 


24 46.238 


- 40 


41 31.128 


15.70 


2008-05-28 


1361 


2 


60 


0.90 


vol 



Note. — Description of columns: (1) Object name; (2) Right ascension {J2000); (3) Declination (J2000); (4) R apparent magnitude 
from catalog; (5) Date of observation; (6) Expos ure time (s); (7) Num ber of spect ra collected; (8) Mean signal-to-noise ratio; (9) Mea n 
seeing during observation (10) References: G04: IGiommi etalj (1200411 ; S94: iStickel et all U994I) : VOl: IVeron-Cetty Veronl j200lD : 
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are reported in Figure [T] while confirmed non- 
BL Lac spectra are shown in Figure [2] All 
BL Lac spectra can be accessed at the website 
http://www.oapd.inaf.it/zbllac/ In particular, 20 
sources have no intrinsic lines, 3 of them show 
weak absorption features arising from their host 
galaxy and other 3 show weak narrow emission 
lines. For only one object (OX 183) a weak broad 
emission line is apparent. We identify spectral 
features arising from the source or from interven- 
ing systems with the atomic species (if applicable) 
while the Diffused Interstellar Bands are marked 
with the label DIB. Telluric absorption are in- 
dicated by the symbol ©. The identification of 
the lines, their equivalent width, full width at half 
maximum and redshift estimation are given in Ta- 
ble [2] for BL Lac objects and in Table [3] for the 
other sources. We also report in Appendix A short 
comments on interesting individual sources . 

3.1. The continuum emission 

In the optical region the emission of BL Lac 
objects is characterised by a superposition of a 
non-thermal component arising by the accreting 
nucleus and a thermal one that originates from its 
host galaxy. In particular, the emission from the 
nucleus can be described by a power law model 
{F\ oc A~") while the galaxy contribution can be 
modeled by an elliptical galaxy te mplate spectrum 



with M]^°** = -22.90 ± 0.50 (ISbarufatti et al 



l2005al ). For the objects that do not present host 
galaxy absorption features we fitted the contin- 
uum adopting a single power law, while for the 
BL Lac where galactic absorption are detected we 
performed the fit by adding the host galaxy tem- 
plate contribution, following procedure described 
in Paper II and III. Resulting power law indexes 
are reported in Table [2] The mean value is 
cxmean = 0.87±0.35 which is consistent with result 
for power law indexes obtained in Paper I, II and 

III (a„ean -0.90 ±0.40). 

3.2. Line detection and EWmin 

We measured the minimum detectable equiva- 
lent width (EW„ii„) for each spectrum adopting 
the recipe described in Paper I. We considered as 
line canditates all the spectral features revealed 
by an automatic procedure selecting EW above 
EW min- Each candidate has been carefully in- 
spected for the consequent feature identification 



or rejection. Spectral parameters, line properties 
and relative identification and EWmi„ for the BL 
Lac objects investigated in this paper are reported 
in Tabled 

3.3. Redshift lower limits 

Despite high quality of the available instru- 
mentation in terms of throughput and the high 
signal-to-noise ratio reached, 20 objects still show 
a featureless continuum. For these sources, we 
applied the procedure fully illustrated in Paper 
II and III in order to measure a lower limit on 
the redshift of each source. Briefly, assuming that 
the BL Lac host galaxy is a giant elliptical with 
j^host ^ _22.90±0.50 if one measures the EW„„ 
and the nucleus apparent magnitude it is possibile 
to infer a lower limit on the z of the source using 
Equation (Paper III) 



EWob 



(1 -I- z) X EWq 
1 + P/Aiz) 



(1) 



where EWohs is the observed minimum equivalent 
width, EWo is the equivalent widt h of the fea- 



ture i n the host galaxy template of iKinnev et al 



p is the nucleus-to-host flux ratio, z is the 
redshift and A{z) is the aperture correction. The 
lower limits are reported in Table [2l 

4. The mean spectrum of BL Lac objects 

In order to investigate the overall spectral char- 
acteristics of the BL Lac objects, we constructed 
a mean optical spectrum combining the normal- 
ized spectra of the 23 sources with z measured 
through intrinsic emission or absorption lines. In 
details, for each spectrum we first carefully re- 
moved spectral region contaminated by telluric, 
interstellar medium or intervening system absorp- 
tion and then we corrected the spectrum to rest 
frame reference by adopting the measured z. We 
finally averaged the spectra according to their S /N 
ratio. The obtained mean spectrum is reported 
in Figure [3] and it covers the spectral range 1800 
A- 6700 A, rest fr ame. The first panel shows 
the mean spectrum assuming for the continuum 
a power-law with index a = 0.90 while the nor- 
malised spectrum is reported in the bottom panel. 
We measure the spectral features detectable on the 
average spectrum reporting the line identification, 
equivalent width and full width at half maximum 
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Table 2 

Spectral line parameters for the BL Lac objects 



Object 


z 


a 


R 


EW„.i„ 


Line ID 


A 




Type FWHM 


EW (observed) 










A 




A 




km 


A 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) 


(11) 


GC 0109+224 


> 0.10 


0.85 


14.20 


0.11 












PKS 0140-059 


> 0.46 


0.45 


18.50 


0.45 












RBS 0337 


0.232 


1.30 


16.00 


0.35 






















Ca II 


4850 


0.233 


g 


-1-0.60 ± 0.15 












Ca II 


4889 


0.232 


g 


-1-0.55 ± 0.15 












G Band 


5304 


0.232 


g 


-1-0.80 ± 0.13 












Mg I 


6375 


0.233 


g 


-1-1.00 ± 0.18 


RBS 0499 


> 0.45 


1.00 


18.60 


0.65 












PKS 0439-299 


> 0.68 


0.30 


20.10 


1.00 






















Mg II (?) 


4708 


[0.68] 


i 


-1-2.15 ± 0.50 


RBS 0589 


> 0.45 


1.20 


17.00 


0.20 












ZS 0506-1-056 


> 0.15 


0.92 


15.30 


0.18 












MS 0622.5-5256 


0.513 


1.00 


19.30 


0.65 






















Ca II 


5953 


0.513 


g 


-1-1.50 ± 0.20 












Ca II 


6002 


0.513 


g 


-1-1.45 ± 0.20 












G Band 


6501 


0.512 


g 


-1-1.40 ± 0.40 


PKS 0735-F17 


> 0.18 


0.84 


15.70 


0.20 












OJ-131 


0.539 


0.44 


17.20 


0.25 






















[O II] 


5736 


0.539 


c 1200 ± 150 


-1.15 ± 0.30 












H /3 


7481 


0.539 


e 1100 ± 100 


-0.70 ± 0.15 












[O III] 


7703 


0.538 


e 1250 ± 100 


-2.50 ±0.20 


RX J1022. 7-0112 


> 0.30 


1.34 


16.60 


0.20 










RX J1023.9-4336 


> 0.24 


1.28 


15.80 


0.12 












PKS B1056-113 


> 0.27 


0.66 


16.50 


0.20 












lES 1106-1-244 


> 0.37 


0.78 


18.20 


0.65 












RX J1117.0-I-2014 


0.140 


1.40 


16.10 


0.40 






















[O II] 


4245 


0.139 


c 1200 ± 120 


-1.20 ± 0.20 












Ca II 


4489 


0.140 


g 


+0.70 ± 0.15 












Ca II 


4522 


0.140 


g 


+0.60 ± 0.10 












G Band 


4908 


0.140 


g 


+0.70 + 0.20 












Na I Blend 


6710 


0.139 


g 


+ 1.10 + 0.40 


RX J1241. 8-1455 


> 0.32 


1.05 


16.80 


0.20 












lES 1239-1-069 


> 0.47 


1.00 


18.10 


0.35 












4C-12.46 


> 0.37 


0.43 


17.60 


0.30 












OQ-240 


> 0.10 


1.00 


14.20 


0.15 












RBS 1478 


> 0.50 


0.95 


18.60 


0.50 












OW-080 


> 0.47 


0.42 


18.10 


0.25 












OX-183 


0.870 


0.16 


18.70 


0.50 






















Mg 11 


5236 


0.870 


e 4400 ± 300 


-5.80 + 1.60 












[O II] 


6970 


0.870 


e 1100 ± 200 


-3.50 + 0.70 


RBS 1899 


> 0.53 


0.93 


18.60 


0.45 










RX J2319. 6-1-1611 


> 0.97 


0.96 


17.20 


0.25 






















M.£! II (?) 


5527 


[0.97' 


i 


J-1,50 + 0.30 


11£S ■^■',22-109 


, II.Kj 


O.Sl 

















Note. — Description of columns: (1) Object name; (2) average redshift (or lower limit); (3) power-law spectral index; (4) R magnitude 
measured on the spectrum; (5) Minimum equivalent width detectable; (6) Line identification (7) Barycenter of the line (8) Redshift of 
the line; (9) Line type e: emission line from BL Lac; g: absorption line from BL Lac host galaxy; i: absorption line from intervening 
system; (10) Full Width at half maximum; (11) Measured equivalent width 
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Table 3 

Spectral line parameters for the objects confirmed non-BL Lacs 



Object 


z 


R 


Line ID 


A 

A 




Type 


FWHM 
km 


EW (observed) 

A 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


PKS 0302-623 


1.350 


17.70 




















C III] 


4485 


1.349 


e 


5500 ± 300 


—20.00 ± 2.00 








Mg II (?) 


5756 


[1.056] 


i 




-1-2.50 ± 0.30 








Mg II 


6585 


1.350 


e 


7000 ± 300 


— 35.00 ± 8.00 


PKS 0410-519 


1.254 


17.10 




















C III] 


4300 


1.252 


e 


5200 ± 300 


-63.50 ± 5.00 








Mg II 


6310 


1.253 


e 


3400 ± 300 


—48.00 ± 6.00 


PKS 0443-387 


0.537 


16.80 




















Mg II 


4298 


0.535 


e 


7400 ± 500 


— 53.50 ± 3.00 








[Ne V] 


5143 


0.537 


e 


1100 ± 100 


—0.60 ± 0.20 








[Ne V] 


5268 


0.537 


e 


1500 ± 200 


-2.60 ± 0.30 








[O II] 


5730 


0.537 


e 


1200 ± 100 


-2.50 ±0.20 








[Ne III] 


5946 


0.537 


e 


2000 ± 300 


-5.30 ± 0.60 








H e 


6098 


0.537 


e 


1200 ± 200 


-1.70 ± 0.30 








H S 


6304 


0.537 


e 


1050 ± 150 


-0.90 ± 0.20 








H 7 


6670 


0.537 


e 


5300 ± 400 


-30.00 ± 3.00 








H/3 


7455 


0.534 


e 


5300 ± 400 


-30.00 ± 3.00 








[O III] 


7697 


0.537 


e 


1100 ± 100 


-35.00 ± 8.00 


PMN J0529-3555 


0.323 


19.20 


















Mg II 


3705 


0.323 


e 


6400 ± 400 


-90.30 ± 10 








Ca II 


5206 


0.323 


g 




-1-4.60 ± 0.70 








Ca II 


5250 


0.323 


g 




-1-4.50 ± 0.50 








G Band 


5689 


0.323 


g 




-1-4.50 ± 0.50 








[O III] 


6560 


0.323 


e 


1300 ± 100 


-4.70 ±0.80 








[O III] 


6623 


0.323 


e 


1100 ± 200 


-9.80 ± 1.20 


PKS 0539-543 


1.191 


17.60 


















C II] 


5097 


1.191 


e 


3200 ± 300 


-1.40 ±0.30 








Mg II 


6135 


1.191 


e 


4000 ± 400 


-50.00 ±8.00 


PKS 0607-605 


1.100 


17.30 




















C III] 


4000 


1.095 


e 


3400 ± 300 


-30.90 ± 3.00 








Mg II 


5880 


1.100 


e 


4000 ± 300 


-45.50 ±4.00 








[O II] 

Mg II (?) 


7821 


1.100 


e 


1100 ± 100 


-13.50 ± 2.00 


OR 103 


1.833 


15.70 


3778 


[0.349] 


i 




+2.80 ± 0.5 








C IV 


4393 


1.834 


e 


5200 ± 300 


-13.50 ± 2.00 








C III] 


5410 


1.833 


e 


4800 ± 400 


-5.00 ±0.70 


PKS 1815-553 


1.633 


20.80 




















C III] 


5030 


1.634 


e 


6800 ± 600 


-92.00 ± 13.00 








Mg II 


7370 


1.632 


e 


6100 ± 400 


-115.00 ± 15.00 



Note. — Description of columns: (1) Object name; (2) average rcdshift; (3) R magnitude measured on the spectrum; (4) Line identi- 
fication (5) Barycenter of the Hne (6) Redshift of the line; (7) Line type e: emission line from BL Lac; g: absorption line from BL Lac 
host galaxy; i: absorption line from intervening system; (8) Full Width at half maximum; (9) Measured equivalent width 
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2000 3000 4000 5000 6000 

Wavelength (A) rest frame 



Fig. 3. — Mean spectrum of BL Lac objects obtained combining the 23 objects of our campaign in which 
intrinsic spectral features are detected. The first panel reports the mean spectrum assuming for the contin- 
uum a power-law with index a — 0.90 (which corresponds to the mean spectral index of the whole BL Lac 
sample). In the second panel normalised spectrum is shown. 
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Table 4 

Spectral line parameters from the BL Lac mean optical spectrum. 



Lino identification A EW FWHM EW [min; max] FWHM [min; max] N„bs / Ntot Line typo 

A A km A km s^^ 

(1) (2) (3) (4) (5) (6) (7) (8) 



C III] 


1909 


1.20 ± 0.10 


4000 ± 350 


[1.30; 1.90] 


[2000; 4000] 


3/4 


o 


Mg II 


2800 


1.00 ± 0.20 


3500 ± 200 


[0.30; 3.20] 


[1900; 6100] 


9/13 


o 


[O II] 


3727 


0.30 ± 0.10 


1300 ± 200 


[0.10; 2.50] 


[1000; 4700] 


9/18 


e 


Ca II 


3934 


0.80 ± 0.20 




[0.30; 6.10] 




12/17 


g 


Ca II 


3968 


0.65 ± 0.10 




[0.40; 4.20] 




12/17 


g 


G Band 


4305 


0.75 ± 0.10 




[0.42; 4.00] 




12/17 


g 


H /3 


4862 


0.25 ± 0.10 




[0.50; 3.00] 




2/14 


g 


[O III] 


5007 


0.50 ± 0.10 


1100 ± 200 


[0.80; 1.47] 


[500; 2400] 


3/14 


e 


Mg I 


5162 


0.90 ± 0.20 




[0.26; 14.00] 




3/12 


g 


Na I 


5892 


0.40 ± 0.10 




[1.00; 3.00] 




2/6 


g 


H a 


6563 


0.90 ± 0.15 


1300 ± 200 




* 


1/3 


e 



Note. — Description of columns: (1) Ion identification; (2) Rest frame wavelength; (3) Measured equivalent width (4) Full width at 
half maximum (5) Equivalent width measurement range (6) Full width at half maximum measurement range (7) Ratio between number 
of detections and number of spectra in which the line could be measures (8) Type of line: e: emission line from BL Lac; g: absorption 
line from BL Lac host galaxy; 



in Table m We note that due to the large difference 
in BL Lacs of the ratio between the non-thermal 
emission and the thermal one associated to stel- 
lar component (their host galaxy) the mean BL 
Lac spectrum cannot be considered as represen- 
tative as that of quasars. On the other hand, it 
illustrates emission or absorption lines that could 
be measured in individual objects. In particular, 
absorption lines from the host galaxy are always 
present in the spectrum but their detectability de- 
pends both on S/N ratio and on the nucleus-to- 
host ratio of the source (see also 15. 4p . Conversely, 
emission lines could be missing. 

5. Discussion and conclusions 
5.1. Summary on the campaign 

According to the 11th Veron-Cetty catalog of 
AGN objects (Veron 2003) at the time of the 
scheduhng of our campaign in 2003, 797 sources 
were known as candidate or confirmed BL Lacs 
and 313 objects of the list were left without a 
reliable measure of their redshift. Considering 
only the sources observable from the southern 
sky {S < 15) a total of 169 BL Lacs objects re- 
mains of unknown redshift. From this set, we 
selected objects with V < 21 in order to se- 
cure an acceptable S/N ratio spectrum for each 
source. This additional constraint allowed us to 



select ~ 155 observable sources from the Paranal 
with VLT-^FORS. We obtained time for the ob- 
servation of 88 sources selected combining the 155 
objects from the 11th Veron catalog with fur- 
ther 44 sources from Padovani and Giommi 1995a. 
Thanks to the high quality optical spectra, we con- 
firmed as bona-fide BL Lacs 69 sources while 15 
have been classified as QSO or FSRQ. The re- 
maining 4 objects are stars. For the sources that 
have been classified as bona-fide BL Lac we were 
able to measured a spectroscopic redshift for 23 
of them. In particular, in 12 sources we clearly 
detected the absorption features from their host 
galaxy and in 15 objects emission lines are ob- 
servable. Broad emission lines (such as Mg II or 
C III]) are detected in only 10 BL Lacs. We mea- 
sured the luminosity of Mg II emission lines in 
all the sources in which the line is detected. The 
mean value is logL(Mg II) ^ 42.40 ± 0.30, where 
L (Mg II) is in erg s~^ . This value is lower of 
an order of magnitude of that of quasars which is 
logL(Mg II) - 43.60 ± 0.45 (IShen et all (l201ll) l 
but the dispersions are rather larger. For the re- 
maining 46 sources we inferred a redshift lower 
limit (see Section 13. 3p . The full dataset, with 
spectra and line identification, is available at the 
website http:/ /www. oapd.inaJ.it/zhllac/. 
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5.2. Optical beaming factor 



In the radio emission, the relativistic beaming 
is manifest in BL Lac objects because of the fre- 
quent appearance of jets with superluminal veloc- 
ities. This is required in high energy emissions (in 
7-rays) by the short time scale of variability, which 
can b e accommodated in a stand ard SSC model 
(e.g. ( Maraschi fc Tavecchiol 2"003() ) only invoking 
large Doppler beaming factors. Following the uni- 
fied model of Active Galatic Nuclei we assume that 
the optical emission arising from the nuclear com- 
ponent is basically the superposition of a thermal 
emission of luminosity Ltu associated with a disk 
structure and a non thermal one Lnth, which dom- 
inates the output power in BL Lac objects. Since 
the former component is isotropic and the latter 
supposedly is beamed, we can de fi ne as beaming 



factor (see e.g. iFarina et al.l (|2012l ). lLandoni et al 

(tml)) 

5=-^ (2) 

Lth 

where Ltot = Lth + Lnth- 

We further assume that the thermal compo- 
nents of BL Lac objects can be described as that 
of quasars. The large collection of Type I AGN 
spectra available in the SDSS allows to find rela- 
tions between the broad line luminosity Lune and 
the contiguous (thermal) luminosity for Type I 
quasars. They are generally expressed in the form 



lof 



XL\{Xref) 



= K 



(3) 



For simplicity, we refer only to Mg II relation since 
it i s the tightest one (K = 1.91±0.26 , A^ef = 3000, 
see lDecarh eral] (|201l[ ): [Shaw et al.l (|2012| )V Tak- 
ing into account the contribution of the cosmolog- 
ical redshift and combining equation (2) and (3) 
with the usual defintion of the equivalent width , 
one has 



d={l + z)x 



A 



ref 



10^' X EW 



(4) 



d is obviously equal to '^1, within the dispersion, 
in the case Type I quasar. When a broad Mg II 
line is observed in a BL Lac object Eq. (4) gives 
an estimate of the beaming. If no emission line is 
detected one can still use Eq. (4) and from the up- 
per limit on the equivalent width deduce a lower 
hmit on S. We suppose that the redshift is such 




Fig. 4. — Optical beaming factor for the sources 
observed in the campaign. Red dashed line in- 
dicates the lower limit value calculated through 
Equation 4 assuming z ~ 0.70. Filled blue circle 
are the 5 of bona-fidc BL Lac objects with Mg 
II emission line while the open magenta circle are 
the beaming factor for other intermediate sources 
(classified as confirmed non-BL Lacs) . Dotted line 
delimit the area of the EW —6 plane where the 
source are of intermediate nature between pure 
QSO {S -1) and BL Lacs (^ > 4). 
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that the relation can be applied because the Mg 
II emission line is within band of observation and 
we also consider the positive 1-a dispersion on the 
relation, which yields K = 2.17 . In particular, for 
featureless BL Lacs we assume z ~ 0.7 while for 
the BL Lac in which the redshift is known we use 
the actual value for the redshift. The wide-spread 
criterion of defining a BL Lac a Type I AGN with 
line EW < 5, which was adopted also in this pa- 
per, therefore translates through Eq. (4) in S > 4. 
In Figure |4] we show the lower limit relation ex- 
pressed in Eq. (4) for the beaming factor of BL 
Lacs with no broad emission lines (dashed line) 
and the actual value of S for the sources in which 
the broad Mg II emission line is detected (filled 
blue circle). We note that ^ 7 BL Lac candidates 
considered in our sample have spectral properties 
that do not completely satisfy the criterion for BL 
Lac classification. On the other hand, these ob- 
jects differ also from the typical spectrum of Type 
I AGN. Interesting enough they have their optical 
beaming factor (open magenta circle of Figure |4]) 
much larger of that of quasar, where 6 should be ~ 
1, suggesting a substantial non-thermal emission 
from the nucleus which indicates that the objects 
are indeed of intermediate nature between BL Lac 
and quasars. 

5.3. Mg II intervening absorption systems 

Since BL Lac objects are extragalactic sources 
it is reasonable to expect that their spectra could 
present absorption featu res from interven in g sys - 



tems. In particular, IStocke fc Rectou ([1997|) 



claimed that a strong excess of a factor ~ 4 on the 
number of Mg II intervening system is apparent 
in the spectra of BL Lac objects, with respect to 
those of quasars. 

In our sample there are 7 confirmed BL Lacs ob- 
jects in which absorption systems are detected. 
These absorptions are interpreted as Mg II inter- 
vening systems on our line of sight toward the BL 
Lac. In particular, the total number of absorb- 
ing systems is Nobs — 8, partly associated with 
BL Lacs with unknown redshift. The most natu- 
ral interpretation is that the features are due to 
intervening Mgll at Aq = 2800. The features asso- 
ciated to intervening Mg II systems are observable 
in our sample in the redshift interval Zmm ~ 0.5 
and Zmax ^ 1.5 . 

The total munbcr of BL Lacs we consider is 



rit = 69 (23 with redshift measured and 46 lower 
limits on z) and those which have z < Zmin are 
13. Therefore, excluding sources with z < Zmin 
the number of BL Lacs which can contribute to 
the observed absorptions is at most 56. On the 
other hand the minimum number which for sure 
contributes to the statistics of the absorbers is 8, 
which are the observed systems associated to BL 
Lacs in our sample, plus the number of BL Lacs 
without intervening absorptions with z > Zmin 



that are 20. FoUowing lZhu fc MenardI (|2012l ). who 
consider QSOs in the SDSS, for an arbitrary point- 
ing direction the normalised mean number of Mgll 
absorbers with EW > 1 A is uq ~ 0.40 in the red- 
shift range given above. Within these assumptions 
we expect that the number of Mg II absorption 
systems is 11< Nexp<22. Since in our sample we 
observe Nobs = 8 interve ning systems we do not 
confirm the suggestion of IStocke fc Recton (|l997t ) 
and we conclude that absorbers related to BL Lacs 
are not more numerous than those associated to 
QSOs. Moreover, if we assume that the distribu- 
tion of Mg II intervening system for BL Lac is not 
different from those of quasars, the number of ob- 
served absorption system could suggest that the 
BL Lacs without z in our sample should have low 
value for their redshift, in order to comply with 
the discussed statistical properties for intervening 
system. 

5.4. Future perspective: BL Lac spec- 
troscopy in the ELT era 

In the spectrum of a BL Lac object the absorp- 
tion lines associated to the stellar components of 
the host galaxy are always present. On the other 
hand, their detectability depends both on the S/N 
ratio and on the nucle us to host ratio (N/H) of 



the source (see e. g. ISbarufatti et all (|2005al) ) 



With the state of the art 8-mt class telescopes 
and modern instrumentation it is possible to ob- 
tain spectra with S/N about 100 and R 2000 for 
object of V ^ 17 (assuming exposure time of ^ 1 
hour). For sources with > 20 the situation is 
even worse since the S/N typically drops below 50. 
This limits the capability to detect the host galaxy 
features either for the sources with high intrinsic 
N/H ratio and for faint objects. For this rea- 
son, in order to significantly improve these limits, 
a future 40-mt class telescope like the European 
Extreme Large Telescope (E-ELT) is mandatory 
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Fig. 5.— First panel: simulation of a BL Lac spectrum with V = 17, S/N = 2000, z = 0.40 and N/H = 
100. Aperture correction is 0.35 . Second panel: simulation for a BL Lac spectrum with V — 20, S/N = 
500, N/H = 40, z — 0.70 and aperture correction of 0.50. For both cases the exposure time is 3600 s and 
the adopted slit is a 1" x 1". The efficiency of the spectrograph is considered similar to the VLT-F0RS2. 
The solid black line is the BL Lac overall spectrum while dashed line is the power law associated to the non 
thermal emission. Blue dotted line is the superposed host galaxy template spectrum. 
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to investigate these kind of objects. 



As example of this capabihty we report here a 
simulation of the expected spectrum of a BL Lac 
object observed by E-ELT. Assuming that the BL 
Lac spectrum is the sum of a non thermal emis- 
sion, modelled with a power law, and a stellar 
component describe d by a giant elliptical galaxy 
(jKinnev et al.lll996[ ) we performed the simulation 



of two cases. We consider an exposure time of 3600 
s and a spectrograph with a dispersion of ~ 1.0 
A I pixel, with an efficiency similar to the VLT- 
F0RS2. We take into account a standard galaxy 
characterised by Mr = —23 and i?e// = 10 kpc. 
We adopted an aperture correction for a slit of 
1" X 1", similarly to the aperture effect correction 
fully described in Paper II. Results are presented 
in Figure [5] . In the first case (top panel) , the 
E-ELT can obtain a spectrum with S/N 2000 
for an object with V = 17 and z = 0.40. The 
features associated to the host galaxy are clearly 
detectable with an intrinsic N/H up to 100 in the 
band 5200 - 6200 A. The second simulation (bot- 
tom panel) is carried out assuming V = 20 and z 
= 0.70. In this case, the S/N ratio of the spectrum 
is ~ 500 and the intrinsic N/H ratio limit reached 
is ^ 40. We note that the use of the E-ELT allows 
to increase the S/N ratio of the spectra of a factor 
~ 20 for the first case and ~ 10 for the latter re- 
spect to typical S/N ratios reached with the VLT. 
In the light of these results, we are confident that 
the advent of 40-mt telescope instrumentation can 
significantly decrease the fraction of BL Lacs that 
are still left without z. 

We acknowledge useful conversations with E. P. 
Farina and A. Veronesi 
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Fig. 1. — Spectra of BLLs sources. Top panel: flux calibrated spectra. Normalized spectra on bottom. 
Telluric bands are indicated by ®, spectral lines are marked by line identification, absorption features from 
interstellar medium of our galaxy are labeled by ISM, diffuse interstellar bands by DIB. The flux density is 
in units of 10~^^ erg cm~^ A~^. 
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Fig. 2. — Spectra of QSOs or confirmed non-BLLs sSSices. Top panel: flux calibrated spectra. Normalized 
spectra on bottom. Telluric bands are indicated by ©, spectral lines are marked by line identification, 
absorption features from interstellar medium of our galaxy are labeled by ISM, diffuse interstellar bands by 
DIB. The flux density is in units of 10~^^ erg cm~^ s~^ A~^. 



PKS 0539-543 



PKS 0607-605 




1.5 



H 1 1 1 1 h 






5000 



6000 



7000 





1.5 




8000 



4000 



5000 



OR 103 




0.3 



0.2 



0.1 




0.8 



4000 5000 6000 

Wavelength (A) 



7000 



2P 

8000 5000 



6000 



7000 80( 
PKS 1815-553 




6000 7000 
Wavelength (A) 



8000 



Fig. 2. — continued 



REFERENCES 

Abdo, A. A., Ackermann, M., Ajcllo, M., et al. 2010, ApJS, 188, 405 
Becker, R. H., White, R. L., & Edwards, A. L. 1991, ApJS, 75, 1 
Berger, K. 2011, International Cosmic Ray Conference, 8, 167 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Costamante, L., & Gliisellini, G. 2002, A&A, 384, 56 
Decarli, R., Dotti, M., & Treves, A. 2011, MNRAS, 413, 39 
Falomo, R., & Ulrich, M.-H. 2000, A&A, 357, 91 

Farina, E. R, Decarli, R., Falomo, R., Treves, A., & Raiteri, C. M. 2012, MNRAS, 424, 393 



Giommi, P., Piranomonte, S., Perri, M., & Padovani, P. 2004, arXiv:astro-ph/0411093 



Giommi, P., Piranomonte, S., Perri, M., & Padovani, P. 2005, A&A, 434, 385 

Healey, S. E., Romani, R. W., Cotter, G., et al. 2008, ApJS, 175, 97 

Kinney, A. L., Calzetti, D., Bohlin, R. C, et al. 1996, ApJ, 467, 38 

Landoni, M., Falomo, R., Treves, A., et al. 2012, A&A, 543, A116 

Majumdar, P. 2011, International Cosmic Ray Conference, 8, 42 

Mao, L. S. 2011, New A, 16, 503 

Maraschi, L., & Tavecchio, F. 2003, ApJ, 593, 667 

Massaro, E., Giommi, P., Leto, C, et al. 2009, A&A, 495, 691 

Nilsson, K., Pasanen, M., Takalo, L. O., et al. 2007, A&A, 475, 199 

Okc, J. B. 1990, AJ, 99, 1621 

Padovani, P., & Giommi, P. 1995, MNRAS, 277, 1477 

Prandini, E., Bonnoli, G., Maraschi, L., Mariotti, M., & Tavecchio, F. 2010, MNRAS, 405, L76 

Rector, T. A., Stocke, J. T., Perlman, E. S., Morris, S. L., & Gioia, I. M. 2000, AJ, 120, 1626 

Rector, T. A., & Stocke, J. T. 2001, AJ, 122, 565 

Sbarufatti, B., Treves, A., & Falomo, R. 2005a, ApJ, 635, 173 

Sbarufatti, B., Treves, A., Falomo, R., et al. 2005b, AJ, 129, 559 Paper I 

Sbarufatti, B., Treves, A., Falomo, R., et al. 2006a, AJ, 132, 1 Paper II 

Sbarufatti, B., Falomo, R., Treves, A., & Kotilainen, J. 2006b, A&A, 457, 35 

Sbarufatti, B., Ciprini, S., Kotilainen, J., et al. 2009, AJ, 137, 337 Paper III 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Seifert, W., Appenzeller, I., Fuertig, W., et al. 2000, Proc. SPIE, 4008, 96 



22 



Shaw, M. S., Romani, R. W., Cotter, G., et al. 2012, ApJ, 748, 49 
Shen, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 45 
Stickel, M., Fried, J. W., & Kuehr, H. 1993, A&AS, 98, 393 
Stickel, M., Meisenheimer, K., & Kuehr, H. 1994, A&AS, 105, 211 
Stocke, J. T., & Rector, T. A. 1997, ApJ, 489, L17 
Veron-Cetty, M.-R, & Veron, P. 1993, A&AS, 100, 521 
Veron, P. 1994, A&A, 283, 802 

Veron-Cetty, M.-P., & Veron, P. 2001, A&A, 374, 92 
Yang, J., & Wang, J. 2010, PASJ, 62, L23 
Zhu, G., & Menard, B. 2012. larXiv:12ir.62T5l 



This 2-column preprint was prepared with the AAS lAT^^X macros v5.2. 



23 



A. Notes on individual sources. 

Adopting the procedure fully described in Pa- 
per II and III we decompose the observed spec- 
trum of objects in which the host galaxy absorp- 
tion features are detected. Briefly the spectral 
decomposition is obtained assuming that the ob- 
served spectrum is a superposition of a non ther- 
mal component, described by a single power law, 
and a thermal emission from a model of a giant 
elliptical galaxy. 



MS 0622.5-5256 is a well known BL Lac 
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Fig. 6. — Spectral decomposition (rest frame) 
for the objects RBS 0337, MS 0622.5-5256, RX 
J1117.0-H2014and RBS 1899. The solid line shows 
the fitted spectrum while dots is the observed one. 



RBS 0337 This BL Lac source has been ob- 



served in the 7-rays band by FERMI (jAbdo et al.l 
2OIOI ) and it has been detected in the MIR/NIR 
bands in th e 2MASS o bservation of BL Lac objects 
II project (iMaol l I2OIII ). Since no imaging studies 
of this object are available, we inferred from the 
spectral decomposition presented in Figure |6] that 
the absolute magnitude of the nucleus and host 
galaxy are respectively M^"*^ 
= -23.80. 



24.30 and Mf 



object (see e .g. iPadovani fc Giommi (Il995l) . 



Rector et al.l (120001)') th at has been studied by 



Sbarufatti et al.1 ( 2005a[) who reports a photomet 
ric redshift of 0.41 through the detection of the 
host galaxy from HST snapshot i mages. Moreover 



further spectroscopical studies (|Sbarufatti et al 



(|2006bh ) report a lower limit of z > 0.49 that 
is obviously consistent with the measurement of 
the redshift obtained in this work. From the 
spectral decomposition obtained (see Figure |6]) 
we measure that the absolute magnitude of the 
nucleus and host galaxy are respectively M^"'^ 
= -23.30 and Mf = -22.70 which are com- 
parable to those meas ured from HST imaging 
Sbarufatti erall (|2005al) . 



OJ 131 This BL Lac object has been studied 
in det ail in the optical band bv iFalomo fc Ulr"icb 
(|2000l) who suggests a lower limit of z > 0.5 from 
the non detection of the host galaxy, assumed to 
be a standard luminous elliptical (A/j^ = —23.80). 
The spectrum of the source presented in this paper 
does not show detectable spectral features deriving 
from the host galaxy. Our new redshift estimate 
is carried out through weak emission lines arising 
from the nucleus itself. 

RX J1117.0+2014 This is a BL Lac source 
studied in a l arge spectral rang e from 7-rays 
with FERMI (lAbdo et all I2OIOI) to the radio 
band ([Becker et al.lll99l[ ). In 2001, this BL Lac 
has been proposed as a TeV candidate source 
(jCostamante fc Ghisellinil 120021 ) but up to now 
no VHE observation are available for this source. 



Moreover, .Nilsson et al.l (|2007l ) measured the ap- 
parent magnitude of the nucleus and its host 
galaxy through imaging decomposition obtaining 



that 



,Host 



— 16.31 ± 0.31. We decomposed the 



observed spectrum (see Figure [6]) in order to de- 
rive the absolute magnitude of the nucleus and its 
host galaxy that are respectively M^""^ — -22.70 
and Mf = -21.50. 

OQ-l-240 is a bright BL Lac object belonging 
to Padovani and Giommi catalog that has been 
observed in 7— rays fro m GeV up to TeV band by 
MAGIC and VERITAS (|Bergeij (jionh . lMaiumdar 
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(|201ll )V Re cent study of thi s sourc e m 7— rays car- 
ried out by IPrandini et aD (|2010l ). based on the 
absorption on Very High Energy photons by the 
Extragalactic Background Light, suggested a red- 
shift of 0.24±0.05 that is consistent with the lower 
limit derived in thi s paper. Moreover, a dopting a 
similar procedure, Yang fc Wang ( 20101 ) propose 
an upper limit for z of 1.19 for this object, well 
consistent with previous results. 



OX 183 is a well kwown BL Lac object (see 
e.g. iRector h Stockd (j200ll) l that has been stud- 
ied through Hubble Spa ce Telescope imaging by 
Sbarufatti et al.l ()2005al ) who suggested a lower 
limit to its redshift of 0.76 that it is obviously con- 
sistent with the z derived from Mg H A5236 and 
[O H] A6970 emission lines reported in this work. 



RBS 1899 Tentatively, there seem to be a Ca H 
break AA 6322 6377 (EW -1.40±0.30, -1.00±0.3) 
in the region around ~ 6300 A that is contami- 
nated by atmospheric absorption, yielding a z = 
0.607. Our spectral decomposition (see Figure H]) 
suggests that M^"'= = -24.60 and Mf = -23.40. 
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